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Abstract 

We present a random walk approximation to fractional Brownian mo- 
tion where the increments of the fractional random walk are defined as a 
weighted sum of the past increments of a Bernoulli random walk. 
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The purpose of this brief note is to describe a discrete approximation to frac- 
tional Brownian motion. The approximation works for all Hurst indices H , but 
take slightly different forms for H < \ and H > ~. There are already several 
discrete approximations to fractional Brownian motion in the literature (see, 
e-g-, HI], [I], 0, HD1, H, 0, 0, [E] for this and related topics), and the ad- 
vantage of the present approach is that the increments of the fractional random 
walk is given as a weighted sum of past increments of an ordinary (Bernoulli) 
random walk. This gives an excellent understanding of the dynamics of the 
process and is a good starting point for stochastic calculus with respect to frac- 
tional Brownian motion. A similar idea is exploited in much greater generality 
by Konstantopoulos and Sakhanenko in [5] , but they assume that H > ^ , while 
the present paper is mainly of interest when H < ^ . 

The discrete approximation is based on Mandelbrot and Van Ness' [5] moving 
frame representation of fractional Brownian motion: 

x t = c H [ ((t - r) H -i - (-r)J-*) db r 
where the scaling constant ch is given by 



,ii ^R i „ A 2 1 A 2 \/T(2H + 1) sin(nH) 

ch=[ (1 + u) 5 - u H -i ) du-> ] - v 







2HJ r(JT + i) 

(see also [3]). This representation will be used to establish the convergence 
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1 The main theorem 



To state the main result, we need some notation. For each natural number N, 
let AtN = jf and think of 

T N = {kAt N | k e 1} 

as a timeline. We let Tjy denote the nonnegative part of T. It is convenient to 
use the following convention for sums over elements in T/v: 

t 

E / W = /( s ) + f( s + At) + • • • + /(t - At N ) 

r—s 

Note that the lower limit s is included in the sum, but the upper limit t is not. 
We shall also write Af(t) — f(t + Aijv) — fit") for the forward increment of / 
at t. 

For all t £ Tjv, let Wjv(i) be independent random variables taking values ±1 
with probability \. We shall write AB N (t) — v / ATjvwjv(i) and think of Bjy as 
a Bernoulli random walk approximating Brownian motion. For < H < 1 and 
N € N, define a process X h ,n : n N x T+ -> R by X/v^(0) = and 

AX h ,jv(s) = 2<:kA^ 5 ABjv(s) + (H--){s-r) H -?At N AB N {r) 

r— — oo 

(using, e.g., Kolmogorov's one series theorem, see [T^], one easily checks that 
the sum converges a.s.) where the constant Kh is defined by 

f -(#-±)C(l-#) for if < i 
[ 1 for H > \ 

(as usual, C(s) = YlnLi n ~ s when s > 1). Except for the Mandelbrot- Van Ness 
scaling factor cjy, Xh,n will be our random walk approximation to fractional 
Brownian motion. For convergence puposes it will be convenient to think of 
Xh.n a s a cadlag process defined on [0, oo), and we do this simply by assuming 
that Xh,n is constant between points in Tjv- 

Remark: Note that the increment AXj^ j jv(s) is a weighted sum of increments 
of the Bernoulli random walk Bn — it is a linear combination of the current 
coin toss wjv(s) and all previous coin tosses ojn{t) 1 r < s. Observe also that 
since lim^i —(H — |)C(| ~ H) = lim s ^i(s — 1)C( S ) = 1) the two cases meet 
continuously at H — i. For H > i we may actually choose Kh as we please 
since the term will vanish in the limit (see below), but Kh — 1 is the natural 
value and probably the one that gives best results in numerical work. 

We are now ready to state the main result. Note that when H = 5, 
AXi N (t) = ABjv(t) and the theorem just reduces to the classical convergence 
of a Bernoulli random walk to Brownian motion. 



2 



Theorem 1 (Main Theorem) For all real numbers H, < H < 1, the pro- 
cesses ChXh,n converge weakly in D([0, oo)) to fractional Brownian motion with 
Hurst index H. 

Notation: In the rest of the paper, we drop the notational dependence on N 
and H, and write simply X, B, T, At for Xh,n, Bn, TV, Atjy etc. when no 
confusion can arise. 

As we are interested in understanding the dynamics of fractional Brownian 
motion, we have defined X by specifying its increments AX(s). To prove the 
main theorem, we need an expression for X(t). This is just a small calculation: 

t t t s * 

X(t) = J2^X(s) = J2K H At H -^AB s + J2 E (H--)(s-r) H -iAtAB r 

s—0 s—0 s—0 r=—oo 

Changing the order of summation, we have 

t t 

X^^KnA^-^Bt + Y, E (H--)(s-r) H -lAtAB r 

r=0 s=r+At 

+ E jZ(H-\){s-r) H -iAtAB r 

r— — oo s—0 

where B t = X)*=o ^ s a ran dom walk converging to Brownian motion. Ob- 
serve that when H > |, the first term Kh At H ~ 2 Bt vanishes when N — > 00 
(this is why the choice of Kh is irrelevant in this case), but when H < |, the 
term explodes. In this case we have a delicate balance between two terms going 
to infinity, and a correct choice of Kh is crucial. 

The idea is now to simplify the expression for X by replacing the sums 
^2(H — |)(s — r) H ^i At by the corresponding integrals J(H — |)(s — r) H ~i ds, 
and then performing the integration. This works nicely for H > 5, but when 
H < i, one of the integrals diverges, and we have to be more careful. Put 
crudely, it is the divergence of this integral that will cancel the divergence of 
the term K H At H -^B t . 

We are ready to prove the main theorem, and start with the simplest case. 

2 The case H > \ 

We start from the expression 

t t 1 
X(t) = At H -iB t +Y / E (H--)(s-r) H -^AtAB r + 

r=0 s=r+At 
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+ E Y,(H-l)(s-r) H -iAtAB r 

r— — oo s— 

above (remember that Kh = 1 in this case). Since H > i, we have no problem 
with convergence, and if we let e^{r,t) be the error term: 

e N (r,t):= £ (H k(s r) H ^ At f (H - ]-)( S - r) H ~Us, 



s=r+A* 

we get 
t 



E (H-h(s-r) H --*At= [ (H-h(.s-r) H -l d S + e N (r,t) = 



s=r+At 



{t-r) H -^ - At H -^ +e N (r,t) 



Similarly, with 

5 N (r, t) := J^iH ~ b( s ~ r) H ~*At - f {H - \){s - r) H ^ ds, 

— n Z Jo ^ 



s=0 

we get 

t 



E( ff - b( s - rf-^At = [ (H- h(s - rf-i ds + S N (r, t) = 
= (t-r) H -^ -(-r) H -^ +S N (r,t) 

This means that 

t 

x(t) = E ((* - r ^ H ~ h + ew w) ASr+ 



+ E ((*-r) H-| -(-r) K " 4 + <Mr)) AB r = 

r— — oo 

t t 

= E ((«-0' i4 -H+ 4 )ABr + ^e W (r,t)AB r + £ 5 N (r,t)AB r 

r— — oo r— r=— oo 

We want to prove that X converges weakly to fractional Brownian motion. 
According to Theorem 1 in 5j, it suffices to show that E(c 2 H X(t) 2 ) — ► t 2H . 
This follows immediately from the Mandelbrot- Van Ness representation and 
the following lemma. 

Lemma 2 For \ < H < 1 : 

fi) E((El= tN(r,t)AB r f) < {H-\)HAt™^ 
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(n) E ((E°=-oo M', t)AB r )i) <(H- I) 2 C(3 - 2H)At 
Proof: (i) We first observe that 

r* 

r+At 



2 if 



e*(r,t) = 2 (H-b(s-r) H -iAt- f (H r) H -i ds > 



s=r+At 



since X)l=r+At(^ — h)( s ~ r ) H * ^ ^ s an u PP er Riemann sum for the integral. 
Since E«=r+2At ~ l)( s ~ r) H ~i At is a lower Riemann sum, we also have 

< e N {r,t) <(H- ^)At H -^At = (H - ^)At H -^ 

Thus ^ t 

E((J2e N (r,t)AB r )A ^Y,e N (r,t) 2 At< 

\ r=0 / r=0 

< J2(H - -) 2 Ai 2H " 1 Ai <(H- -) 2 tAt 2H - x 
(ii) Using approximating Riemann sums as in part (i), we see that 
0<5 N (r,t)<(H-±)(-r) H -?At, 

and thus 

e(( J2 S N (r,t)AB r A=j2S N (r,t) 2 At< ]T (H \f{-rf H ^ At 3 

\ r— — oo / — oo r— — oo 

Letting r = —kAt, we get 

£ Mr,t)AB r ) 2> ) <f>-i) 2 fc 2ff - 3 At 2ff = 

\ r= — oo / fc=0 

= (if - i) 2 C(3-2 J ff)At 2H 

This completes the proof of the lemma (and also the proof of the Main Theorem 
for the case H > |). □ 

3 The case H < \ 

Again we start from the expression 

t t 

X^^KnA^-iBt+Y, Yl (H - -)(s - rf-i AtAB r + 

r=0 s=r+At 
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r— — oo s—0 



In this case, one of the integrals we worked with above diverges, and we have 
to be more careful. Let us start with a closer look at the term ^2 s =r+At(^ — 
|)(s — r) H ~iAt. We obviously have 



oo _. OO _. 

E (H - -)(s r) H ~iAt £(ff - 2 )(S r)H ~ lAt 

s=r+At ' s=t 

and if we let r = NAt, s = kAt, we get 

oo OO 

E (H--)(s-r) H -iAt= E (^-2)(fcAt-7VAi) ff -iAt 

s=r+Ai fc=7V+l 

1 oo 1 oo 

= (tf-±)Af*-* E (k-N) H ~i = (H-±)At H -iJ2 nH ~ i 

k=N+l n=l 

= (H- l)At H -?((l -H) = -K H At H -? 
Substituting this into the expression for X(t), we get 

*(*) = £ £ ~( H - b( s r) H -^AtAB r 



r=0 s=t 
t 



+ EE(tf-^)(*-r)"- | AiAB r 



r— oo s—0 



The two sums in this expression have less dangerous limits than the one we just 
got rid of, and can be approximated by integrals. If we let 



1 P OO I 

r(r, t) := E ~(H - ^ § Ai - / " o^ S - ds ' 

we get (remember that < 5 ): 

°° -| /'OO 1 

£-(#- 2)(«-^ _i Ai= / -{H--)(s-r) H -U S + l N {r,t) 

s=t J* 

1 -i s—oc- 1 

-(s-r) H -i +e N (r,t) = (t - r) H -? + e N (r, t) 
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Similarly, if wc let 



~S N (r, t) := J2 ~(H \)(s - rf^At f -{H - rf-i ds, 

s=0 "'° 



we get 



E(# - |)( s - r )^ § At = /V - ^)( s - r ) H ~ % ds - *) 



= (a-r)""* " _^( r ,t) = (t-r)*-*-(-r) ff -*-for(r,f) 
L J s=0 

We thus have 

t 

X W =£((*" r ) H_ " + *)) A£? r + 



+ E ((* - r) H ~^ - (—r) H ~i — S N (r, t)J AB r 

= E ((t-r)*"* -(-r);-*)**+ 



r= — oo 

i 



E ^( r ' t)AB r — S N (r,t)AB r 



r— — oo 



To prove that c# converges weakly to fractional Brownian motion, we can 
now longer use Theorem 1 of [5] as in the previous case since this theorem 
requires that H > ~. However, the first term in the expression above obviously 
converges weakly to 



((t-r)*"i-(-r);-*) 



and the next lemma shows that error terms go uniformly to zero. Using the 
Mandelbrot- Van Ness representation, we then get the Main Theorem for H < |. 

Lemma 3 For each H , < H < |, there is a constant Kh € (independent 
of N and t ) such that 



X(t)- E ((t-r) H ---{-rf + -- 



AB r 



< K H At H 



Proof: It clearly suffices to show that there are constants Ch,Dh € K+ (inde- 
pendent of N and t) such that 



r=0 



< C* H Ai ff and 



S N (r)AB 7 



r— — oo 
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We begin with the e^-case. By definition 
e N (r,t) = J2-(H--)(s-r) H -iAt- / -{H - -){s - r) H ~Us 

s=t Jt 

Since YlT=t ~ (H — ^)( s — r) H ~i At is an upper Riemann sum for the integral 
f t °° -(H - |)( fl - r)*"* da, and £r= t+ A t -(# - - r)*-*Ai is a lower 
Riemann sum, we have 

0<e>(r,i)<-(^-^)(t-r) H -fAi 

Hence (remember that |A£? r | = At 2) 



J2^n{t, t)AB, 



r=0 



<j2-(H-b(t-r) H -iAti 



r=0 



If we let t = KAt, r = kAt, we can rewrite the last sum as 
V -(H \){K k) H ~iAt H < —{H \)C{\ - 



fc=0 



This completes the ejy-part of the argument. 

Turning to the term Y^t=-oo ^n{t)AB t , we first observe that by definition 

~Sn(v) ^-{H- h{s rf-iAt f —(H his r) H -i ds 
s =o 2 Jo 2 

Again, X)*=o — (H ~ l)i s ~ r) H ~i At is an upper Riemann sum, and we easily 
see that 

< S N (r) < —{H - h(-r) H ~^At 



Letting r = —kAt, we get 



E 



S N (r)AB r 



^ E -{H-\)(-r) H -iAti< 



1 00 1 ■? 

< -{H I) At" £ = -{H \)C{1 H)At H 



fc=0 



This proves the lemma (and hence the Main Theorem for the remaining case 



H<\). 



□ 
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